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THEORETICAL  BUCKLING  LOADS  OF  BORON/ALUMINUM  AND  GRAPHITE/RESIN 
FIBER-COMPOSITE  ANISOTROPIC  PLATES 
by  Christos  C.  Chamis 
Lewis  Research  Center 

SUMMARY 

Theoretical  results  are  presented  for  designing  with  boron/aluminum  composites 
when  these  composites  are  buckling  critical.  The  composites  are  assumed  to  be  rec- 
tai^ular  plates  with  four  simply  supported  edges.  They  are  subjected  to  single  and  com¬ 
bined  normal  and  shear  loads  in  the  plane  of  the  plate.  The  plates  are  made  from  a 
unidirectional  composite  whose  fiber  direction  is  oriented  at  some  angle  to  the  load 
direction. 

The  design  data  are  presented  in  nondimensional  form  as  buckling  load  against 
orientation  angle  for  several  plate  aspect  ratios.  The  results  indicate  that  the  buckling 
loads  of  boron/aluminum  plates  are  independent  of  fiber  orientation  if  the  plate  aspect 
ratio  is  greater  than  approximately  1.  The  bucklii^  loads  are  moderately  dependent  on 
the  orientation  angle  for  plates  with  aspect  ratios  less  than  about  1.  The  buckling  load 
is  independent  of  aspect  ratio  in  plates  with  aspect  ratios  greater  than  about  2. 

Comparison  of  buckling  results  for  boron/aluminum  composite  plates  and 
Thor nel-7 5/epoxy  composite  plates  indicates  that  the  boron/aluminum  composite  plates 
resist  bucklii^  loads  more  efficiently  than  the  Thornel-75/epoxy  composite  on  the  basis 
of  specific  buckling  strength.  The  results  also  indicate  that  the  bucklii^  loads  of  boron/ 
aluminum  composite  plates  can  be  predicted  using  orthotropic  theory  if  their  aspect  ra¬ 
tio  is  greater  than  about  1. 

The  numerical  algorithm  used  to  solve  the  buckling  problem  and  listing  of  the  cor¬ 
responding  computer  code  through  which  the  results  were  obtained  are  included. 


INTRODUCTION 

Feasibility  studies  for  the  space  shuttle  indicate  that  the  use  of  advanced  fiber  com¬ 
posite  structural  components  can  result  in  a  considerable  increase  in  payload  in  the 
shuttle  system.  Boron/aluminum  and  graphite/resin  fiber  composites  are  leading  con¬ 
tenders  for  shuttle  applications  because  these  composites  offer  h^h  stiffness-to-density 


and  high  strength-to-density  ratios.  Panels  made  from  these  materials  will  have  to  meet 
both  material  strength  and  buckling  requirements.  This  report  deals  with  a  theoretical 
investigation  of  the  buckling  of  flat  rectangular  panels  made  from  boron/aluminum  and 
Thornel"7 5/epoxy  fiber  composites. 

Several  papers  have  appeared  recently  dealing  with  the  buckling  of  anisotropic  plates 
(refs.  1  to  8).  However,  design  data  are  not  available  for  flat  panels  subjected  to  com¬ 
pressive  loads  and  made  from  advanced  fiber  composites  such  as  boron/aluminum  and 
graphite/resin.  The  method  described  in  reference  3  has  proved  efficient  in  buckling 
studies  of  boron/epoxy  plates.  This  method  is  used  herein  to  generate  design  data  for 
boron/aluminum  plates.  Some  data  for  Thornel-75  graphite/epoxy  resin  plates  are  also 
generated  for  comparison  purposes.  Data  were  generated  for  aspect  ratios  of  1/2,  1, 

2,  and  4. 

The  panels  considered  are  anisotropic  and  simply  supported.  They  are  subjected  to 
combined  in- plane  (normal  and  shear)  load  (fig.  1).  The  material  is  a  unidirectional 
composite  with  the  fiber  direction  oriented  at  an  arbitrary  angle  to  the  load  direction 
(fig.  1).  The  analytical  algorithm  used  is  the  assumed  mode  method  in  conjunction  with 
the  Galerkin  method.  A  computer  code  was  developed  based  on  the  Galerkin  method,  and 
the  code  was  used  to  generate  the  theoretical  design  data  presented  herein.  A  brief  de¬ 
scription  of  the  analytical  method  is  given  in  the  report.  All  symbols  are  defined  in  ap¬ 
pendix  A.  The  numerical  algorithm  used  to  solve  the  resulting  eigenvalue  problem  is 
described  in  appendix  B.  Input  data  sample  sheets  with  explanations  are  given  in  appen¬ 
dix  C.  A  listing  of  the  computer  program  with  sample  cases  is  given  in  appendix  D. 


y 


Figure  1.  -  Fiber-composite  plate  geometry  an^  loading  -  all  four  edges  simply  supported 
(aspect  ratio  =  alb). 
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BRIEF  DESCRIPTION  OF  UNDERLYING  THEORY 


The  underlyii^  theory  for  buckling  loads  of  anisotropic  plates  is  described  in  refer¬ 
ence  3  with  pertinent  discussions  in  references  6  to  8.  Briefly,  this  theory  consists  of 
expressing  the  potential  energy  of  a  plate  in  terms  of  displacement  variables.  Taking 
the  variation  of  the  potential  energy  function  yields  the  field  equation  and  the  correspond¬ 
ing  boundary  conditions.  The  resulting  system  then  is  solved  by  the  assumed  mode 
technique  in  conjvmction  with  the  Galerkin  method. 

The  equation  resultii^  after  the  variation  of  the  energy  function  is 
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(The  notation  is  defined  in  appendix  A. )  The  area  integral  represents  the  field  equation, 
and  the  line  integrals  represent  the  boundary  conditions. 

The  assumed  buckling  mode  described  in  reference  3  is  represented  by  a  Fourier 
double  sine  series.  This  mode  satisfies  the  imposed  boimdary  conditions,  but  it  does 
not  satisfy  the  natural  boundary  conditions  if  the  material  and  structural  axes  do  not  co¬ 
incide.  However,  the  mode  is  forced  to  satisfy  the  natural  boundary  conditions  approx¬ 
imately  through  the  Galerkin  method  as  discussed  in  reference  3. 

Substituting  the  assumed  mode  in  equation  (1),  applying  the  Galerkin  method,  and 
carrying  out  the  algebra  result  in  a  set  of  linear  equations  which  represent  the  eigen¬ 
value  problem  of  the  plate.  This  system  is  coupled  for  either  a  combination  of  shear 
and  normal  loads  and/or  noncoincident  material  and  structural  axes. 

The  eigenvalue  problem  is  solved  by  using  the  Power  method,  which  is  a  highly  ef¬ 
fective  iterative  numerical  technique  in  seekir^  the  largest  eigenvalue  of  the  system. 
The  indicial  equations  which  were  used  to  generate  this  system  and  the  Power  method 
are  given  in  appendix  B  in  outline  form. 
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BRIEF  DESCRIPTION  OF  COMPUTER  PROGRAM 


The  numerical  algorithm  described  in  appendix  B  has  been  transformed  into  a  com¬ 
puter  code  which  is  rather  simple  and  can  be  generated  from  the  information  supplied 
in  appendix  B.  A  FORTRAN  IV  compiled  listing  of  a  computer  program  based  on  the 
algorithm  in  appendix  B  is  given  in  appendix  D,  with  sample  cases  and  output.  Input 
data  sample  sheets  are  given  in  appendix  C. 

The  inputs  to  the  code  are  composite  system  identification,  fiber  volume  ratio, 
orientation  angle,  plate  aspect  ratio,  and  flexural  rigidities.  The  outputs  are  the  num¬ 
ber  of  terms  in  the  assumed  mode  series  expansion  required  for  convergence,  the  rela¬ 
tive  error  between  the  last  two  iteration  cycles,  the  bucklii^  load,  and  topo-plot  data 
of  the  buckled  shape  of  the  plate  normalized  with  respect  to  the  largest  deflection. 

The  algorithm  described  in  appendix  B  runs  into  difficulty  when  the  shear  buckling 
load  of  a  plate  is  sought.  In  this  particular  case,  the  difficulty  is  bypassed  by  includii^ 
normal  loads  which  are  a  very  small  fraction  of  the  shear  load.  Further  discussion  on 
why  these  difficulties  arise  is  presented  in  reference  3. 

THEORETICAL  DESIGN  DATA 

The  theoretical  design  data  generated  herein  are  based  on  the  schematic  illustrated 
in  figure  1.  In  this  figure,  the  type  of  loadii^  condition,  the  plate  geometry,  and  the 
fiber  orientation  are  defined.  The  x-y  coordinate  reference  system  is  referred  to  as 
the  structural  axes  system.  The  fiber  direction  coordinate  system  which  is  located  at 
the  angle  d  from  the  structural  axes  system  is  referred  to  as  the  material  axes  sys¬ 
tem.  The  loading  conditions  are  identified  by  Ny,  and  N^y  as  is  noted  in  the 
figure. 

The  flexural  stiffnesses  required  in  calculating  the  buckling  loads  were  calculated 
by  usii^  the  computer  code  described  in  reference  9.  Typical  values  of  the  elastic  con¬ 
stants  of  the  plate  along  its  material  axes  are  given  in  table  I  for  boron/aluminum  and 
Thornel- 75/epoxy  composites  with  a  fiber  volume  ratio  of  0.  5.  The  flexural  rigidities 
are  computed  as  functions  of  the  orientation  angle  using  the  data  in  table  I. 
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TABLE  I.  -  THEORETICAL  UNIDIRECTIONAL  COMPOSITE 


PROPERTIES  AT  FIBER  VOLUME  RATIO  OF  0.  5 


Property 

Boron/aluminum 

Thornel-75/  epoxy 

Longitudinal  modulus, 

24.  2X10®  (35.  OxlO®) 

26.  OXIO®  (37. 8X10®) 

N/cm^  (psi) 

6.9X10®  (I.OXIO®) 

Transverse  modulus, 

16.  8X10®  (24.  3X10®) 

N/cm^  (psi) 

8.0x10®  (11.6x10®) 

0.  44x10®  (0.  63x10®) 

Shear  modulus, 

N/cm^  (psi) 

Poisson’s  ratio 

0.  24 

0.25 

Density,  g/cm 

2.62  (0.095) 

1.  55  (0.  056) 

(Ib/in. 

Buckling  of  Boron/Aluminum  and  Thornel -75/Epoxy  Panels 

Buckling  loads  for  a  single  loading  condition  for  panels  made  of  boron/aluminum  and 
Thornel-75/epoxy  are  illustrated  in  figure  2,  where  the  specific  buckling  stress  has  been 
plotted  as  a  function  of  the  orientation  angle  for  an  aspect  ratio  of  2.  The  schematic  in 
the  figure  illustrates  the  type  of  load  condition  as  well  as  the  orientation  angle.  As  can 
be  seen,  in  this  figure,  boron/aluminum  composites  are  more  efficiently  utilized  than 


Figure  2.  -  Specific  buckling  stress  of  two  fiber-composite  plates  - 
all  four  edges  simply  supported.  Fiber  volume  ratio,  0.5;  aspect 
ratio  a/b,  2. 
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Figure  3.  ^  Buckling  loads  for  two  fiber-composite  plates  -  all  four  edges  simply 
supported.  Fiber  volume  ratio,  0.5;  aspe-t  ratio  a/b,  2. 


Thornel-75/resin  composites  in  structures  which  are  critical  in  bucklir^  as  measured  by 
the  specific  buckling  stress. 

Buckling  load  comparisons  where  the  panels  are  loaded  in  the  y-direction  are  il- 
lustrated  in  figure  3.  In  this  figure,  it  can  be  seen  that  the  boron/aluminum  composite 
panel  is  considerably  stronger  in  buckling  than  the  corresponding  Thornel-75/epoxy 
panel.  In  this  plot,  the  nondimensional  buckling  load  parameter  is  plotted  as  a  function 
of  the  orientation  angle  for  the  fixed  panel  aspect  ratio  of  2. 

Buckling  Loads  for  Individual  Loading  Conditions 

Design  data  for  boron/aluminum  panels  which  are  subjected  to  compressive  load  in 
the  x-direction  are  illustrated  in  figure  4  as  a  function  of  orientation  for  various  aspect 
ratios.  The  important  point  to  be  noted  from  this  figure  is  that  the  buckling  load  is  in¬ 
dependent,  or  almost  independent,  of  orientation  angle  in  panels  where  the  aspect  ratio 
is  approximately  greater  than  1.  Another  point  to  be  noted  is  that  the  buckling  load  de¬ 
pends  only  moderately  on  the  orientation  angle  in  panels  of  aspect  ratio  less  than  1. 

A  cross-plot  of  figure  4  is  illustrated  in  figure  5.  The  nondimensional  load  is 
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Figure  4.  -  Buckling  loads  for  boron/aluminum  composite  plates,  with  all 
four  edges  simply  supported,  subjected  to  normal  (Nj^)  compressive  load. 
Fiber  volume  ratio,  0.5. 


ORIENTATION 

ANGLE, 


Figure  5.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected  to 
normal  compressive  load.  Fiber  volume  ratio,  0.5. 
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plotted  as  a  function  of  panel  aspect  ratio  for  various  orientation  angles.  The  insensi¬ 
tivity  of  the  buckling  load  as  a  function  of  orientation  angle  in  panels  with  aspect  ratios 
greater  than  about  1  is  clearly  illustrated  in  this  figure.  The  dashed  line  symbol  is  used 
to  represent  these  curves  to  emphasize  that  buckling  load  values  were  computed  only  at 
the  aspect  ratios  1/2,  1,  2,  and  4. 

Buckling  loads  for  panels  which  are  loaded  in  the  y-direction  only  are  given  in  fig¬ 
ure  6.  As  can  be  seen  in  this  figure,  the  buckling  loads  are  almost  independent  of  the 
orientation  angle  in  panels  with  aspect  ratios  of  1/2  and  greater.  The  buckling  load,  on 
the  other  hand,  is  very  sensitive  to  the  aspect  ratio  in  panels  with  aspect  ratios  of  ap¬ 
proximately  2  or  less. 

Buckling  loads  for  a  boron/alxxminum  composite  panel  loaded  with  shear  only  are  il¬ 
lustrated  in  figure  7.  The  points  to  be  noted  in  this  figure  are  the  following: 

(1)  There  is  a  mild  buckling  load  dependence  on  the  orientation  angle  for  panel 
aspect  ratios  of  less  than  about  1. 

(2)  The  buckling  load  is  relatively  independent  of  orientation  angle  for  panel  aspect 
ratios  of  greater  than  about  1. 

(3)  The  buckling  load  is  very  sensitive  to  the  aspect  ratio  in  panels  with  aspect  ra¬ 
tios  less  than  2,  and  this  dependence  becomes  rather  insignificant  for  panel  aspect  ra¬ 
tios  greater  than  2. 


ASPEa 


Figure  6.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  ail  four  edges  simply  supported,  subjected 
to  normal  (Ny)  compressive  load.  Fiber  volume  ratio,  0.5. 


Figure  7.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected 
to  shear  load.  Fiber  volume  ratio,  0.5. 

Ay 
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Buckling  Loads  for  Two  Equal  Simultaneous  Loadings 


Buckling  loads,  when  the  panel  is  loaded  with  equal  loads  in  the  x-  and  y- 
directions,  are  shown  in  figure  8.  In  this  figure,  the  nondimensional  buckling  load  pa¬ 
rameter  is  plotted  as  a  function  of  orientation  angle  for  various  panel  aspect  ratios. 
The  results  in  this  figure  show  that  the  buckling  load  is  slightly  dependent  on  the 


Figure  8.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected 
to  combined  normal  compressive  (Nj^  =  Ny)  loads. 

Fiber  volume  ratio,  0.5.  ^ 


orientation  angle  for  panel  aspect  ratios  less  than  1,  and  practically  independent  of 
orientation  angle  for  aspect  ratios  greater  than  1.  The  buckling  load  is  sensitive  to  panel 
aspect  ratio  for  aspect  ratios  less  than  or  equal  to  2.  This  dependence  becomes  insig¬ 
nificant  for  panel  aspect  ratios  greater  than  2.  The  curves  of  the  buckling  load  as  a 
function  of  the  independent  variables  indicated  in  figure  8  parallel  the  curves  of  the 
buckling  loads  indicated  in  figures  4  and  6,  for  the  individual  loadings. 

Buckling  loads  for  the  case  when  the  panel  is  loaded  in  the  x-direction  combined 
with  shear  are  shown  in  figure  9.  The  curves  of  the  buckling  load  for  this  loading  con¬ 
dition  are  parallel  to  those  of  the  individual  cases  (figs.  4  and  7).  Buckling  loads  for 
the  case  when  the  panel  is  loaded  in  the  y-direction  combined  with  shear  are  illustrated 
in  figure  10.  The  buckling  load  in  this  figure  seems  to  be  practically  independent  of 
orientation  angle  for  the  aspect  ratios  investigated.  However,  it  is  quite  sensitive  to 
the  panel  aspect  ratio  for  aspect  ratios  less  than  2. 
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NONDIMENSIONAL  BUCKLING  LOAD,  N^rb  /h  NONDIMENSIONAL  BUCKLING  LOAD, 


Buckling  Loads  for  Three  Equal  Simultaneous  Loadings 


Buckling  loads  for  panels  which  are  loaded  in  the  x-  and  y-directions  combined  with 
shear  are  shown  in  figure  11.  The  schematic  in  this  figure  indicates  the  type  of  loadings 
and  their  respective  ratios.  The  nondimensional  buckling  load  is  plotted  as  a  function  of 
orientation  angle  for  various  panel  aspect  ratios. 

Comparing  corresponding  curves  from  figures  8  and  11,  it  is  seen  that  the  addition 


y 


Figure  11.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected  to 
combined  normal  (Ny  =  and  shear  (N^^y  =  loads. 
Fiber  volume  ratio,  ^0.5.  ^ 


of  the  shearing  load  decreases  the  buckling  load  of  the  panel  only  slightly.  The  point  to 
be  noted  then  is  that  a  panel  subjected  to  compressive  loads  in  the  x-  and  y-directions 
will  resist  almost  an  equal  amount  of  shearing  load  for  approximately  the  same  buckling 
load. 


Buckling  Loads  for  Two  or  Three  Unequal  Simultaneous  Loadings 

Buckling  loads  for  panels  which  are  subjected  to  unequal  loads  in  the  x-  and  y- 
directions  are  shown  in  figure  12.  The  type  of  loading  condition  and  respective  loading 
m^nitudes  are  illustrated  in  the  sketch  given  in  the  figure.  In  this  figure,  the  nondi¬ 
mensional  buckling  load  parameter  is  plotted  as  a  function  of  the  orientation  angle  for 
various  aspect  ratios.  The  curves  of  the  buckling  load  for  this  type  of  loading  condition 


11 


Figure  12.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected  to 
combined  normal  (Nw  =  (l/2)Ny)  loads.  Fiber  volume  ratio, 

n  c  y 


ASPECT 


Figure  13.  -  Buckling  loads  for  boron/aluminum  composite 
plates,  with  all  four  edges  simply  supported,  subjected  to 
combined  normal  (Ny  =  (1/2)NJ  and  shear  (Nj^  =  (l/2)Nj^) 
loads.  Fiber  volumeTatio,  O.b. 

parallel  those  of  the  cases  with  equal  loading  condition,  as  was  previously  discussed  (see 
fig.  8).  One  additional  point  to  be  noted  is  that  the  buckling  loads  of  panels  with  aspect 
ratios  of  greater  than  about  2  remain  almost  invariant  as  a  function  of  aspect  ratio  when 
the  orientation  angle  is  greater  than  about  45°. 

Buckling  data  for  panels  loaded  with  unequal  combined  loadii^  conditions  are  illus¬ 
trated  in  figure  13.  The  loading  condition  for  the  panel  is  illustrated  in  the  schematic 
in  the  figure.  The  buckling  load  is  plotted  as  a  function  of  orientation  angle  for  various 
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aspect  ratios.  A  point  to  note  in  this  figure  is  that,  at  some  orientation  angles,  panels 
with  aspect  ratios  greater  than  2  could  have  greater  buckling  loads  than  panels  with 
aspect  ratios  of  2.  As  can  be  seen,  the  buckling  load  for  a  panel  with  an  aspect  ratio  of 
approxin>ately  2  is  lower  than  for  the  panel  with  an  aspect  ratio  of  4  when  the  orientation 
angle  is  approximately  less  than  45° 

All  the  buckling  data  presented  and  discussed  indicate  that  the  buckling  load  is  insen¬ 
sitive  to  orientation  angle  for  panels  with  high  aspect  ratios.  The  buckling  load  is  mildly 
sensitive  to  the  orientation  angle  in  panels  with  low  aspect  ratios.  This  observation 
leads  to  the  important  conclusion  that  the  buckling  loads  of  boron/aluminum  composite 
anisotropic  plates  can  be  approximately  determined  by  using  classical  orthotropic  theory. 
This  conclusion  is  indeed  a  useful  one,  since  the  buckling  of  orthotropic  plates  has  been 
extensively  treated  in  reference  5. 


POSSIBLE  EXTRAPOUTIONS  OF  DESIGN  DATA 

The  design  data  presented  and  discussed  were  based  on  a  fixed  fiber  volume  ratio  of 
0.  5.  The  data  presented  herein  can  be  used  to  extrapolate  buckling  loads  for  plates  made 
from  composites  with  different  fiber  volume  ratios.  The  results  will  be  very  close  if  the 
variation  of  the  fiber  volume  ratio  is  within  approximately  ±20  percent  of  the  0.  5  value 
which  was  used  in  generating  the  design  data. 

This  extrapolation  is  justified  since  the  buckling  load  is  nondimensionalized  with 
respect  to  both  composite  longitudinal  modulus  and  thickness.  It  is  well  known  that  both 
composite  modulus  and  thickness  depend  on  the  fiber  volume  ratio  (ref.  10),  and  that  this 
dependence  is  approximately  linear  in  the  fiber-volume-ratio  range  0.  4  to  0.  6.  In  this 
sense,  then,  the  extrapolation  using  the  design  data  presented  herein  for  fiber  volume 
ratios  within  ±20  percent  of  0.  5  should  yield  reasonable  results. 

The  computer  code  appended  in  appendix  C  can  be  slightly  modified  to  compute  the 
first  natural  frequency  of  anisotropic  boron/aluminum  composite  panels.  In  reference  5, 
the  vibration  problem  is  discussed,  and  the  analogy  in  computing  the  buckling  load 
and  the  natural  frequency  is  made. 


CONCLUSIONS 

The  discussion  of  the  theoretical  design  data  presented  leads  to  the  following  con¬ 
clusions: 

1.  Design  data  for  the  buckling  of  unidirectional  boron/aliiminum  panels  with  fibers 
oriented  at  any  ai^le  to  the  load  direction  have  been  generated  and  are  reported  herein. 
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2.  Specific  buckling  stress  comparisons  showed  that,  in  general,  boron/aluminum 
composite  panels  are  more  efficient  than  high-modulus  graphite/resin  composite  panels. 

3.  The  buckling  load  of  boron/aluminum  unidirectional  panels  is  practically  inde¬ 
pendent  of  fiber  direction  at  high  aspect  ratio  values.  At  these  aspect  ratios,  the  plate 
can  be  assumed  to  have  its  material  axes  coincide  with  its  structural  axes.  Conse¬ 
quently,  the  classical  buckling  theory  of  orthotropic  plates  can  be  used  to  predict  the 
buckling  load. 

4.  The  buckling  loads  of  boron/aluminum  unidirectional  panels  are  only  moderately 
dependent  on  fiber  direction  at  plate  aspect  ratios  less  than  1. 

5.  Boron/aluminum  composite  panels  loaded  by  normal  in-plane  loads  which  are 
near  the  critical  load  can  carry  considerable  shear  load  before  they  buckle. 

6.  The  buckling  loads  of  boron/aluminum  panels  are  practically  independent  of 
aspect  ratio  at  aspect  ratios  greater  than  about  3. 

7.  The  buckling  loads  of  panels  with  fiber  volume  ratios  within  approximately 
i20  percent  of  0.  5  can  be  extrapolated  from  the  design  data  presented  herein  by  using 
the  appropriate  panel  thickness  and  the  appropriate  composite  longitudinal  modulus. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  September  3,  1971, 

129-03. 
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APPENDIX  A 


SYMBOLS 


a  panel  dimension,  x-direction 

b  panel  dimension,  y-direction 

D  matrix  of  flexural  rigidities,  eq.  (1) 

h  thickness 

K  plate  stiffness  matrix,  elements  given  by  eqs.  (Bl)  and  (B4) 

L  load  matrix,  elements  given  by  eqs.  (Bl)  and  (B4) 

M  summation  index  limit  on  m 

m  summation  index 

N  summation  index  limit  on  n;  with  subscripts,  applied  load 

N  parameter  in  eq.  (Bl) 

Ncr  buckling  load 

n  summation  index 

p  summation  index 

q  summation  index 

W  vector  of  coefficients  in  the  displacement  mode  expansion,  eq.  (1) 

w  displacement  in  z-direction 

X,  y,  z  panel  structural  axes 

1,  2,  3  panel  material  axes 

6  variation 

convergence  tolerances 

9  orientation  angle  of  material  axes  with  respect  to  structural  axes 

X  buckling  parameter 

p  ratios  defined  by  eq.  (B5),  also  weight  density 

a  stress  and  matrix  of  stresses 

Subscripts: 

cr  critical  or  bucklii^  condition 

i  row  index 


] 


column  index 


I  undirectional  composite  property 

m  buckling  modes,  x-direction,  axes 

n  buckling  modes,  y-direction 

r  iteration  cycle 

X,  y,  z  directions  associated  with  the  respective  structural  axes 
1,  2,  3  directions  associated  with  the  respective  material  axes 
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APPENDIX  B 


DESCRIPTION  OF  NUMERICAL  ALGORITHM 

The  numerical  a^orithm  seeks  the  eigenvalue  of  the  following  matrix  equation; 

[K]{W}  =  N[L]{W}  (Bl) 

where  [K]  and  [L]  are  (M  X  N)  x  (M  X  N)  square  matrices,  {W}  is  a  column  matrix  con¬ 
taining  the  and  N  is  defined  subsequently.  The  (i,  j)  elements  of  the  [K]  and  [L] 

matrices  are  given  by  the  following  indicial  expressions: 
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where 


p 

^xy  N 


(B5) 


where  N  is  the  buckling  parameter.  Here  N^,  Ny,  and  N^y  carry  their  algebraic 
signs,  wherea.s  N  is  always  taken  as  a  positive  quantity. 

The  buckling  loads  are  determined  by  finding  the  largest  eigenvalue  in  equation  (Bl) 
using  the  Power  method.  Several  methods  are  available  for  finding  eigenvalues  of  linear 
systems  (ref.  3).  The  Power  method  is  relatively  easy  to  program  and  is  applicable  to 
nonsymmetric  matrices  which  have  real  eigenvalues.  The  method  yields  the  largest 
eigenvalue  and  the  corresponding  eigenvector.  To  apply  the  Power  method,  equation  (Bl) 
is  expressed  in  the  following  form: 


X{W}  =[K]“^[L]{W}  (B6) 

where  X  =  1/N.  The  solution  is  obtained  in  an  iterative  fashion  as  follows: 


(B7) 


where  {W}  is  normalized  relative  to  its  largest  element.  Thus 


where 


{W} 

{W}  = - 

A  \ 


X^  =  max.  (W^)  (i  =  1  (1)  MX  N) 


(B8) 


(B9) 


The  iteration  process  stops  when 


^  1  -  K 

r+1  r 


<€ 


r+1 


(BIO) 


where  e’  is  usually  taken  10'^  <  e’  <  10"®.  This  procedure  converges  rapidly  as  long 
as  the  next  largest  eigenvalue  is  not  close  to  X.  The  Power  method  runs  into  difficul- 
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ties  when  the  shear  buckling  load  is  sought  because  for  this  case  the  magnitude  of  the 
two  largest  eigenvalues  is  the  same.  This  difficulty  is  easily  overcome  by  using  very 
small  values  for  one  or  both  of  the  normal  loads  in  combination  with  the  shear  load  as 
is  described  in  reference  3. 

The  buckling  load  is  obtained  from  the  relation 


1^ 

N 


X 


r+1 


and 


(BID 


Thus  the  smallest  buckling  load  is  obtained  since  is  the  largest  eigenvalue  of 

equation  (B6). 

The  buckling  load  of  the  plate  then  is  determined  by  incrementing  M  and  N  in 
equation  (Bl)  and  applying  the  Power  method  to  compute  N  until  the  convergence 
criterion 


N(M  +  1)(N  +  1)  (B12) 

N(MXN) 

is  satisfied.  The  parameter  e”  is  usually  chosen  as  an  acceptable  percentage  of  N. 
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APPENDIX  C 


INPUT  DATA 
Explanations 


Explanations  appear  in  parentheses  leadered  to  correspondii^  card. 


1  2  i  U 

b  7  8  9  10  II  IZ 
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mm 

■B 
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mm 
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mm 
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Hi 
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^H 
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■mi 
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0. 

0. 
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p 
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0. 
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'  '-1 . 

0 
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0- 

IHH 
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0.5 

.  A’X 

0 . 0 

o’ . 

' 

I'.O 

0. 5 

'-i'- 

1.0 
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l» 
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. 
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■ 
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mm 

BP 

BS 
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Si 
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HHH 

hhhh 

■Bi 
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BH 

BH 

HHH 

hhh 

Bmi 
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NASA-C-836  (REV.  9-14-59] 
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Sample  Input  Data  for  0. 02-Inch-Thick  Boron/Aluminum  Plate 


1  2  S  6  5 

6 

7  B  9  10  II  12 

IS  M.  15  16  17  IS 

19  20  21  22  23  24 

25  26  27  28  29  50 

31  32  35  34  55  36 

49  50  51  52  53  54 

55  56  57  56  59  60 

61  62  6364  65  66 

67  68  69  70  71  72 

75  74  75  76  77  76  79  8' 
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.5 
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MM 
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1 
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42 
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.5.0 
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mm 
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mm 
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BB 
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APPENDIX  D 


COMPUTER  OUTPUT 


C....  BUCKLING  GF  ANISOTROPIC  PLATES 
CCNKCN  /EIGVEC/  W(81) 

CCRKCN  /aRFILl/  01  3, 3  ):/ ARF I L 2/RHO(  3 )  / ARR AY/X C  81 , 81 ) 

REAL  N1,N2,LAMB2 

OIRE^SIC^  VdOflOVt  TilO),  RST(  20f  3  I  ,  T  ITLE<  3 ) 
data  PI/3. 1415927/ 

MGCE^l 

ReAC(5,23)  TITLE 
READ  (5»12)  NP,MU 
READ  (5,13)  A,B 

READ  (5,14)  ((RSTiI,J):,J=l,3),I  =  l,NP) 

1  READ  (5,25)  ( (0(  UJ  ) ,  J=1 ,3) ,  1^1 ,3  )  ,THETA 
DC  11  II*1,NP 

WRITE  (6,15) 

DC  2  J=l,3 

2  RHC(  J):=RST(II ,  J) 

N1=0*C 

CALL  TI^^LFT  CQl) 

3 

call  ^^FULl  (M,A,B) 

CALL  ^INVA  CMM,KSIG) 

IF  (KSIG.EG.2)  GO  TO  10 
CALL  (^FILL2  (M,  A,B) 

CALL  EIGFIN  (NMUAMB2) 

IF  (KSIG.eC.3)  iFLAG^r 

N2==l,0/L/SMe2 

EPS-ABS((N2-N1)/N2) 

IF  (EPS.LE.l.OE-4)  GO  TO  5 

IF  (^^.EC.WU)  GO  TO  4 

SPS=EPS 

W=R+1 

N1  =  N2 

GO  TC  3 

4  WRITE  {6,16) 

5  IF  (IFLAG.EQ.l)  WRITE  (6,17) 

CALL  Tlf^LFT  (02) 

IC=(Gl-C2)/60. 

WRITE(6,24)  TITLE, THETA 
WRITE  (6,18)  M,SPS,EPS,IO 
V^^AX=0.0 

DC  7  1=1,10 
DC  7  J=l,10 
VU,J1  =  0.0 
DC  6  K=1,M 
DC  6  L=1,M 
IJ=(K-1)=«v  +  l 

6  V(  I  ,  J)=V(  I  ,  J  )+V(  TJ  X'SINIFLOATfK*!  )*PI/10.0)*SIN(  FL0AT(L<^J)»PI/10.  ) 

7  V^AX  =  AMAX1 ( VMAX, ABS( V( I ♦ J)  )  ) 

DC  8  1=1,10 

DC  8  J=l,10 
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V(I,J)=V(I,J)/VMAX 

WRITE  (6,19)  A,R,N2,W0DE,RH0,D 

DC  9  1=1,10 

J=1C-I+1 

T(n=Cf.O 

WRITE  (6,20)  J,T(I),(V(K,J),K=1,10) 

1=0 

XI=C.C 

WRITE  (6,20)  I,XI,(T(J),J=1,10) 

WRITE  (6,21) 

GG  TC  11 
WRITE  (6,22) 

CCNTINUE 
GC  TC  1 

FCRRAT  (212) 

FCRRAT  (2F10.4) 

FCRRAT  (9F8.3) 

FCRRAT  (IHl) 

FCRVAT  (23H  MAXIMUM  INDEX  ATTAINED) 

format  (41H  ILL-CONDITIONED  MATRIX  -  RESULT  IN  DOUBT) 

FORMAT  (5X,2HM=,  I1,‘5X,20HPREVIOUS  REL.  ERROR=  ,  IPE 10 . 3 , 5X  ,  1 6HLAST  R 
lEL.  ERR0R=,1PE10.3,5X,.33H  TIME  REQUIRED  FOR  THIS  CASE  WAS  ,I3,8H  S 
2ECCNCS) 

FORMAT  (///7X,1HA,13X,-1HB,  11X,3HNCR,  11X,4HM|0DE,9X,5HRH0-X,9X,5HRH0 
1-Y,9X,6HRHC-XY//3(2X,1PE10.3,2X ) , 6X , 1 1 , 7X, 3 ( 2X , IPE 10. 3, 2X ) ///37H  T 
2HE  C-ARRfY  COLUMNWISE  IS  AS  FOLLOWS/ /9 ( 2X, IPE 10 . 3, 2 X ) ////50X , 13HBU 
3CKLEC  SHAPE//6X, IHY) 
format  (//‘5X,I2,.11(5X,F5.2)) 

FORMAT  (//15X,1H0,9X, 1H1,9X, 1H2,9X,1H3,9X, 1H4, 9X , 1H5 , 9X , 1H6 , 9X , 1H7 
1,9X,1H8,9X,1H9,8X,2H10,6X, IHX) 

FORMAT  (37H  MATRIX  SINGULAR  -  PROCESS  TERMINATED) 

FCRMAT(3A6) 

FCRMAT(/20H  COMPOSITE  SYSTEM  =  ,3A6,12H  AT  THETA  =  ,F5.1) 

FORMAT  (10F8.3) 

ENC 


SUERCUTINE  MFILLl  (MM,A,B) 

COMMON  /ARFILl/  0(3,3) 

COMMON  /ARRAY/  X(81,81) 

INTEGER  F,C 

DATA  PI/3.  141  5927/ 

DC  3  M=1,MM 

DC  3  N=1  ,MM 

DC  3  P=1,MM 

DC  3  C=1,MM 

I=(M-l)*MM+N 

J=  (P-1) ♦MM  +  C 

IF  (  I  .NE.  J )  GO  TO  1 

IF  (P.NE.M.OR.O.NE.N)  GO  TO  2 

X{I,J)  =  Pl»*4*(B*FLnAT(M*>i‘4)*C(l,l)/A**3  +  2.0*FLnAT(M*M*N*N)>!‘(D(l,2) 
1+'2.C#C(  3 ,3  )  )  /  (  A*P  )+A*FLnAT(N**4  )*D(  2, 2)/B«*3)  /4.0 
GC  TC  3 


o  o 


1  TF  ( iwrn  ( ( M  + p )  ^ ( N -f 0 )  •  2  1  •  FQ •  ^  2 

X(  I  .  jr=-8  .0*PI<'*2#D(  U3)’«'FLnAT(  M«N*P*Q^(  2*m«m/(P*P-M*M)  +  1  )/(0’9^0-N« 
IN)  )/A’»’!'2-8. 0*r>I*«2*0(  2,3  )*FLnAT(M*N*P*Q*(2*N't'N/(0*Q-N*N  )+!)/{  P*P-M 

2*K) ) /B**2 
GC  TC  3 

2  X(I,J)=0.0 

3  CCNTIMJE 

return 

ENC 


SUBRCUTINE  I^INVA  (M,KSIG) 

C . MATRIX  INVERSION  BY  GAUSS-JORDAN  ELIMINATION 

CCMMCN  /iRRAY/  X(81,81) 

DIMENSION  Y(81,81),  Z(2,81t81),  K(81) 

equivalence  (Y,Z) 

double  precision  a,r,proo,y 

DATA  CELT,EPS,L0CPS/0.0001,l.0E-8,2/ 

KSIG=C 

IF  (M.NE.l)  GO  TC  3 
PRCi:  =  X(  1,1 ) 

IF  (PROD. NE. 0.0)  GC  TO  2 

1  KSIG=KSIG+2 
RETURN 

2  Xn,l)  =  1.0/PR00 
RETURN 

3  PRrC=1.0 

MM=M-1 

DC  A  1  =  1, M 
K(n  =  i 
DC  A  J=1,M 

4  Y(I,J)=X(1 ,J) 

c . begin  by  FINDING  LARGEST  PIVOTAL  ELEMENT 

DC  9  1  =  1,  M 
A  =  C.O 

DC  5  J=I,M 

IF  ( A8S ( Y ( J , 1 ) ) • LE.  A  )  GO  TO  5 
A=AeS(Y(J, 1)  ) 

L=J 

5  CONTINUE 

IF  (A. EC. 0.0)  GO  TO  1 

c. ...  .rearrange  rows  and  order  array 

N=K(L) 

K(L)=K(  I  ) 

K(  I  )=N 
DC  6  J=1,M 
A=Y(  I,J) 

Y(I,Jl=Y(L,J) 

Y(L  ,  J)  =  A 

REDUCE  PIVCTAL  ROW 
A  =  Y (  I  ,1  ) 

DC  7  J=1,MM 
7  Y( I ,J)=Y(I ,J+1)/A 

Y(  I,M)  =  1.0/fl 
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C . RECtCB  REMAINING  ROWS 

DC  9  L=1,M 
IF  (L.FC.I  )  GO  TO  9 
) 

DC  8  N=1,MM 

Y(LfN)=Y(L,N+l)-A*Y( I,N) 

8  IF  (ABS(V(L»N)).LT.(ABS(Y(L,N+1) )*EPSn  Y(L,NI=0*C 
YU*M)=-^«Y<UM) 

9  CONTINUE 

C . UNSCRAMBLE  INVERTED  MATRIX 

DC  n  1=1, M 

IF  (Kd  ).EC.I  ).  GO  TO  13 
DC  10  J=I,M 

10  IF  <K(J).EC.n  GO  TO  11 

GO  TC  1 

11  DC  12  L=1,M 
A=Y(Ld  ) 

Y(L,  I)  =  Y{Lf J) 

12  Y(L,J1=A 
KU)=K(  I  ) 

13  CONTINUE 

C . OBTAIN  EPRCR  MATRIX 

DO  18  N=l, LOOPS 
TEST=C.O 
DC  15  1=1, M 
DC  15  J=1,M 
R  =  C.O 

DC  14  L=1,M 

14  R=R-7Cl,L,J)^Xn ,L) 

IF  (UEC.JI  R=R+1.0 
TEST=AMAX1(TEST, ABS(RL) 

15  Z{2,UJ)=R 
DC  17  1=1, M 
DC  17  J=1,M 
A  =  0.0 

DC  16  L=1,M 

16  A=A4Z(1,I,L)^Z(2,L,J) 

17  Z(1,I,J)^Z(1, I,J  )+A 

IF  (TEST.LE.DELT  )  GO  TO  19 

18  continue 

KSIG=KS IG+3 

C..*. .TRANSFER  FINAL  INVERSE 

19  DC  20  1=1, M 
DC  20  J=1,M 

20  X(I,J1  =  Z(1,  I,J) 

RETURN 

END 


SLBRCLTINE  MFILL2  (VM,A,B) 
CC^^MCN  /ARFIL2/  RHO(3) 

CCMMCN  /ARRAY/  X (81,81) 

dimension  L(81,81),  C(81,81) 

real  L 

INTEGER  F,C 

DATA  PI /3 . 1415927/ 
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DC  3 

DC  3  N=l,^^ 

DC  3  P=l,^r' 

DG  3  C=ltPP 

I=(P-1»«M^+N 

J= (P-1) «^M+C 

IT  (I.NE.J)  GO  TO  1 

IF  (P.NE.M.OP.O.NE.N)  GO  TO  2 

L(  I,J)  =  -FI**2*(B*FLOAT(M>CM)*RHO(  1) /A+A*FLOAT(N*N  )*RHn(  2  ) /B) /4.0 
GC  TC  3 

1  IF  (PGD(  ((^  +  P)<‘(N+Q),2).EQ.O)  GO  TO  2 

LT I , J )  =  -RM0(3 )  *FLO  AT( 8*M*N*P*Q/ • t  P*P-M*M ) ♦ ( Q*Q-N*N ) ) ) 

GO  TC  3 

2  L<I,J)=0.0 

3  CONTINUE 

DC  5  I  =  1,NN' 

DC  5  J=ltNN 

SUP=C.O 

DC  A  IJ=ltNN 

4  SUP=SUM+X{  ItIJ)*L(  IJ,J) 

5  C(I,J).=  SIP 
DC  6  1=1, NN 
DC  6  J=1,NN 

6  X(I,J)=C(I,J) 

RETURN 

END 


SUBRCCTINE  EIGFIN  (MP,,LAMR2) 
real  LAPei,LAMB2 
CCPNCN  /EIGVEC/  W(81) 

CCPPCN  /ARRAY/  X(8l,8l) 

DIPENSION  V(81) 

LAPPl=O.C 
W( 11=0.5 

IF  (PP.EC.l)  GO  TO  2 
DC  1  1  =  2, PP 

1  W(I)=C.5 

2  DC  A  1  =  1, PP 
SL'P  =  C.O 

DC  3  J=1,PP 

3  SUP  =  SUM  +  X ( I , J1*W ( J  ) 

4  V(I)=SUP 
LAPB2=A6S{V(1)) 

IF  (PP.EC. 1)  GO  TC  6 
DC  5  1=2, PP 

5  LAPe2  =  APAXl(LAMB2, ABSIVI  I)  )  ) 

DC  7  1  =  1, PP 
Vvd  )=V(  I)/LAMB2 

IF  (APS  (  (LAPP.2-LAPR1)/LAMB2)  .LE.l.OE-4)  RETURN 
LAKPlnL APB2 
GC  TC  2 
END 
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CGMPCSITE  SYSTEM  =  BORON/ ALUM  I NUM  AT  THETA  =  30.0 

Nn5  PREVIOUS  REL.  ERRORS  1.302E~03  LAST  REL.  ERROR=  2.519E-03  TIME  REQUIRED  FOR  THIS  CASE  WAS  11  SECONDS 
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CCMPCSITE  SYSTEf^  »  BORON/ ALUM  INUM  AT  THFTA  30,0 

PREVIOUS  REL.  ERROR*  3.292E-02  LAST  REL.  ERROR*  2.30BE-03 
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